Long-range odd triplet superconductivity in SF structures with Neel walls 
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We consider a multidomain superconductor/ferromagnet (SF) structure with an in-plane magne- 
tization, assuming that the neighboring domains are separated by the Neel domain walls. We show 
that an odd triplet long-range component arises in the domain walls and spreads into domains over 
a long distance of the order £t = y D/2nT (in the dirty limit). The density of states variation in 
the domains due to this component changes over distances of the order £t and turns to zero in the 
middle of domains if the magnetization rotates in the same direction in all domain walls. 

PACS numbers: 74.45.+C, 74.78.Fk, 74.50.+r, 75.70.Kw 



I. INTRODUCTION 

The past decade has seen a rapid growth of interest 
in the study of hybrid superconductor - ferr omagnet (SF) 
structures (see, for example, reviews IH l2l3l) . The interest 
in such systems originates from the possibility of finding 
new physical phenomena as well from the hope of con- 
structing new devices based on these structures. New 
physical phenomena arising in these systems are the re- 
sult of a nontrivial interplay of competing types of order- 
ing in superconductors and ferromagnets. Superconduct- 
ing correlations lead in superconductors to the appear- 
ance of the Cooper pairs, that is, the pairs of electrons 
with opposite spins. On the opposite, the exchange in- 
teraction in ferromagnets tries to align the electron spins 
in one direction. In SF structures these two types of in- 
teractions are spatially separated and can coexist despite 
much greater value of the exchange energy h in compar- 
ison to the superconducting order parameter A. 

Due to the proximity effect the superconducting cor- 
relations penetrate the ferromagnet in SF structures. 
The opposite effect, i.e., the penetration of a magnetic 
moment M into the superconductor, also takes place. 
It turns out that the magnetic moment Ms is induced 
in the superconductor. The magnetic moment Mj is 
aligned in the direction opposite to the magnetization di- 
rection of free electrons in the ferromagnet and spreads 
over a distance of the o rder of the superconducting cor- 
relation length £s = y/Ds/A (in the dirty limit) £ On 
the other hand the condensate wave function / pene- 
trates the ferromagnet with an uniform magnetization 
over a much shorter di stance of the order of the 
"exchange length" £ h = Dp jhw& The condensate 
wave function decays in F in a nonmonotonic way as 
f(x) ~ exp(— x/£h) cos(x/£h); it oscillates in space and 
decreases exponentially. This nonmonotonic behavior of 
f(x) leads to a nonmonotonic dependence of the criti- 
cal temperature T c of the superconducting transition in 
SF bilayers and multilayer s 1 ! 2 ! 6 ! 7 ! 8 ! 9 ! 10 and to a 7r-state in 
SFS Josephson junctions, 1 ! 2 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 

In the case of a nonuniform magnetization in the ferro- 



magnet a new phenomenon appears: a triplet component 
of the condensate wave function / (generally speaking the 
condensate wave function is a matrix in the particle-hole 
and spin space) arises in the SF system^ This triplet 
component is an odd function of the Matsubara fre- 
quency to (while the conventional BCS singlet component 
of / is an even function of w) and spreads in the ferromag- 
net over a long distance of the order of £t = \J Df/2ttT. 
This long-range triplet odd-frequency component was 
predicted to exist in a SF structure with nonhomogeneous 
magnetization in Ref. |T7] and this prediction was con- 
firmed for a slightly different case in Ref. In Refs. |n] 
and 0, a SF structure with a domain wall at the SF 
interface was considered, that is, it was assumed that 
Mj? = M(0, sina(a;), cos a{xy\ , where a(x) = Qx in the 
interval < x < ciq and a(x) = Qclq at x > clq. The 
triplet component Jl was shown to arise in the domain 
wall and to penetrate the ferromagnet over a long dis- 
tance £t- Unlike the triplet component in superfluid 3 He 
and in Sr2Ru04, this odd triplet component corresponds 
to s-wave correlations and hence is symmetric in the mo- 
mentum space; therefore it is not destroyed by scattering 
on ordinary, nonmagnetic impurities, and survives in the 
dirty limit. We call this component the long-range triplet 
component (LRTC). The LRTC may also arise in a SF 
structure with a uniform magnetization and spin-active 
interfaced Note that from a macroscopic point of view 
a domain wall at the SF interface can also be consid- 
ered as a "spin-active interface". The LRTC may arise 
in a multilayered SF structure with noncollinear orien- 
tations of the magnetization vector M^ in different Fj 
layers^SiS 1 " In particular, a new type of superconductiv- 
ity (odd triplet superconductivity) has been predicted in 
such structures if the thickness of the F layers d obeys the 
condition: £/j <C d < £t- In this case the Josephson cou- 
pling between neighboring F layers is realized only via 
the LRTC because the singlet component decays very 
fast in the F layers. Therefore superconductivity in the 
transverse direction is due to the LRTC, whereas in-plane 
superconductivity is caused mainly by the singlet BCS 
component. The influence of the LRTC on the critical 
temperature of the superconducting transition in FSF 
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structures with a noncollinear magnetization orientation 
was studied in Ref. |U 

Historically, the odd-frequency triplet pairing was con- 
jectured in 1974 by Berezinsii 23 as a possible mecha- 
nism for superfluidity in 3 He. It turned out later that 
in 3 He another type of triplet pairing (even in frequency, 
odd in momentum) is realized. O dd- freque ncy pairing 
in solids was also studied in Refs. 12 12, '1 2(1 A triplet 
odd-frequency pairing was investigated in Refs. (two- 
dimensional electron gas with repulsion in the presence 
of impurities) andl25l (a Kondo lattice model). A singlet 
odd-frequency pairing was analyzed^ as a possible type 
of pairing in high-temperature superconductors. 

Although several experimental results may be inter- 
preted in terms of the LRTC; 2 Ti 2 ?i 2 ?i 3 ?i 31 i 32 there are 
still no direct experimental evidences in favor of the odd 
triplet superconductivity. Therefore there is a need to in- 
vestigate, in more detail, a possibility to observe this new 
type of superconductivity. One of the important and in- 
herent features of ferromagnets is the domain structure. 
The domain structure in a ferromagnet may essentially 
alter the properties of SF structures. For example, mod- 
ulation of the phase difference tp between two supercon- 
ductors due to internal magnetic fields in the ferromag- 
netic domains, can lead to a negative critical current I c in 
the Josephson SFS junction^ The influence of domains 
on I r in SFS Josephson junctions was studied in Refs. 
and |35[ In Ref. yj the magnetization vector is assumed 
to be in-plane and to rotate around the direction nor- 
mal to the plane of the junction. In this case the LRTC 
arises, and the possible 7r-state may be suppressed due 
to an effective averaging of the exchange field. In Ref. 
a SFS junction with two in-plane domains of opposite 
orientation was studied (no LRTC arises in this case). 
It was shown that if the thicknesses of the domains are 
equal, the critical current I c is always positive. 

In realistic domain structures, the domains are sepa- 
rated by domain walls. Below we shall discuss the case 
when the magnetization vector Mp lies in the plane of 
the F film and the domain walls are of the Neel typei^i 
A limiting case of such a structure was considered in 
Refs. M and |H It was assumed that the vector Mp 
rotates in space continuously 

M F {y) = Mo(0,sinQy,cosQy) (1) 

(we choose the x axis normal to the plane of the F film, 
whereas in Ref. [3^ the z axis is normal to the plane of the 
F film). The possibility of a cryptoferromagnetic state 
in SF structures with the magnetization Mp (y) given by 
Eq. was studied in Ref.|33. The F layer was supposed 
to be very thin: d « ((,. A solution for the Eilenberger 
equation has been found near the critical temperature T c 
of the superconducting transition. It was established that 
in a certain interval of parameters (Q, d, etc.) a homo- 
geneous state in the ferromagnet becomes energetically 
unfavorable and the nonhomogeneous magnetization de- 
termined by Eq. arises in the F film. In Ref. |3^ 
a SF bilayer with the F film of arbitrary thickness was 




FIG. 1: SF systems considered in the paper, (a) SF bilayer 
with half-infinite S and F parts. The domain (y < 0) and 
the region with rotating magnetization (y > 0) are also half- 
infinite, (b) Multidomain F layer of thickness d in contact 
with a bulk superconductor. Depending on the relative ori- 
entation of rotating magnetizations in neighboring domain 
walls, we distinguish the cases of positive and negative chiral- 
ity (Q has the same or opposite sign in the neighboring do- 
main walls, respectively). The proportion between the widths 
of domains and domain walls is chosen only for drawing pur- 
poses. 



studied in the dirty limit. Analysis of a solution for the 
Usadel equation shows that in the case of magnetization, 
uniformly rotating along the F film [see Eq. QJ], the 
condensate function / penetrates into the F film over a 
short distance of the order £?,, . The LRTC is absent in 
this case. According to Refs. lnllllllTl2ll the LRTC ap- 
pears if Mf rotates across the F layer (i.e., depends not 
on y but on x). 

In the present paper we consider a domain structure in 
a thin F film, where domains with antiparallcl in-plane 
magnetizations are separated by the Neel walls (while the 
magnetization does not change across the thin F film). 
This domain structure is realized in real ferromagnetic 
filmsiSi The yz plane is chosen to be parallel to the SF 
interface (see Fig. QJ. We show that the LRTC arises 
at the Neel domain walls and decays exponentially away 
from the domain walls and the SF interface over a long 
distance £r- We calculate the density of states (DOS) 
variation bv in the ferromagnet caused by the proximity 
effect and find that 5v(y) turns to zero in the middle of 
domains in the case of positive chirality. 

The paper is organized as follows. In Sec.^ we formu- 
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late the Usadel equations, the corresponding boundary 
conditions, and investigate the main features of the long- 
range triplet superconducting component which appears 
due to the presence of Neel domain walls. The analysis 
is made on the simplified model with half-infinite S and 
F layers and only one half-infinite region with rotating 
magnetization. In Sec. IIIII we consider the case of the 
multidomain F layer, employing the results of Sec.[H] To 
make the model realistic, in Sec. II VI we take into account 
a finite thickness of the ferromagnet. To illustrate the 
results for the LRTC, we study the density of states due 
to it. Finally, we present our conclusions in Sec. W\ 



II. BASIC EQUATIONS. HALF-INFINITE 
DOMAIN 

Consider a ferromagnet (x > 0) in contact with a su- 
perconductor (x < 0). We assume that the in-plane ex- 
change field h(y) in the F layer, which is proportional to 
the magnetization Mp, depends on y: h(y) — h(Q, 0, 1) 
at y < and h(y) — h(0, sin a(y), cos a(y)) with a{y) = 
Qy at y > 0. This means that the magnetization vector 
Mj? is oriented along the z axis at y < and rotates in 
the yz plane at y > 0. The region with rotating magne- 
tization models a Neel domain wall. The structure that 
we discuss first is shown in Fig. ^a) and contains only 
one half-infinite domain and one half-infinite region with 
rotating magnetization; the thickness of the ferromagnet 
is also infinite. Then we shall use the obtained results to 
describe a realistic structure depicted in Fig.QJb). 

Our goal is to find the condensate Green functions in 
the ferromagnet, induced due to the proximity effect. We 
consider the dirty limit, which means, in particular, that 
hr <C 1, where r is the momentum relaxation time due 
to elastic scattering. 

We employ a widely used model: an exchange field h 
acts on free electrons and there is no attractive interac- 
tion in the F layer leading to the superconducting order 
parameter However, the condensate (Gor'kov) func- 
tions are finite in the F region due to the boundary con- 
ditions at the SF interface. We are interested in distances 
much larger than the Fermi wavelength, therefore we can 
use the quasiclassical Green functions g R ^ A ^> . The matrix 
retarded (advanced) functions the ferromagnet 

obey the Usadel equation^ 



DX7(g"^>X7g 



(2) 

where e is energy, h = h{r^aj, cos a(y) + tq&2 sin a(yj) , 
and a(y) = at y ^ while a(y) = Qy at y ^ 0. 
We assume that the diffusion constants D for electrons 
with spin up and down are equal (this is correct if the 
exchange energy h is much less than the Fermi energy 
Ep). The matrix Green functions g are 4x4 matrices in 
the Gor'kov-Nambu and spin spaces, fj and Oj are the 
Pauli matrices in the Gor'kov-Nambu and spin spaces, 
respectively. 



We represent g R( - A ) in the form 

g R W = ± f 3 a + /, 



(3) 



where the first term is the retarded (advanced) quasiclas- 
sical Green function in the normal state. The supercon- 
ducting correlations, described by f(x,y), are assumed 
to be weak due to the finite interface transparency (re- 
sulting from an oxide barrier or from mismatch in the 
Fermi surfaces of S and F materials). In the considered 
case of weak proximity effect (|/| <C 1), Eq. J2J can be 
linearized. In the Matsubara representation it acquires 
the form 

V 2 / - 2k 2 J - ik\ sgnw({<7 3 , /} cosa(y) 



+ f 3 [<7 2 ,/]sina(y)) = 0, (4) 



where u = ttT(2u + 1), k 2 = \w\/D, k 2 = h/D, the 
square brackets denote the commutator, and the braces 
denote the anticommutator. 

The Green function in the bulk of the superconductor 
is fs — T-2<Jzfs, with fs = A/\/lu 2 + A 2 . We shall use 
the following boundary condition for / at the SF interface 
(x = 0): 



df 

dx 



(5) 



where 7^ = Rb<r, while a is the conductivity of the 
ferromagnet and Rb is the interface resistance per unit 
area. This boundary condition follows from the general 
one o 39 i if two assumptions are made: (1) the proximity 
effect is weak (i.e., 1) and (2) the bulk solution 

fs in the superconductor is unperturbed and valid up to 
the interface (i.e., jb/£s ^> v/<?s)- 

We can rewrite Eq. Q in the form taking into account 
boundary condition (|5|): 

g^+Qyi- 2P J - ik l sgaw({d-3, /} cosa(y) 



+ t 3 [ct 2 ,/] sina(y)) = -— S(x). (6) 

'76 



We must seek for an even solution of this equation, then 
this is equivalent to the problem with the boundary con- 
dition (we have reflected / with respect to x = and now 
solve the equation at all x). However, the requirement 
that the solution is even, will be automatically satisfied: 
as we shall see below, the Fourier harmonics (over x) de- 
pend only on k 2 , hence they are even in k, which means 
that f(x) is even in x. 

Performing the Fourier transformation /(fc, y) 
j dxf(x, y) exp(— ikx), we obtain 

-l-(k 2 + 2fc 2 )/ - ik\ sgnu,({<7 3 , /} cosa(y) 



dy 



+ h[a 2 ,f]sma(y)) = (7) 
'lb 
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At y > the function a(y) is y-dependent, while at 
y < we have a — 0. In the region of positive y one can 
exclude the y-dependence from Eq. with the aid of 
rotation 

f = Uf u U+, (8) 
where U = exp(iT3<7ia(y)/2). As a result, we get (y > 0) 

+ iQf 3 {o-i,^} -ifc 2 sgnw{o- 3 ,/„} = - — (9) 
dy 7 h 

in terms of the new function y). The same equation 
is valid for y < if we set Q = 0: 



<9y 2 



(fc 2 + 2^)/ tl -jfc2 sgnw | ( j 3j / u } 
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(10) 



The original functions / and df jdy are continuous at 
y = 0. Therefore the rotated functions obey the following 
boundary conditions at y = 0: 



M-o) = fu(+Q), 



(11) 



-^y— = -^- + *2 T3 ^ 1 ' / »>- (12) 

Thus we have to solve the linear matrix Eqs. JJjJ (y > 
0) and HI U| l (y < 0) of the second order with the boundary 
conditions and l|12|l at y — 0. We can represent the 
solution in the form 

f u = F(Q)9(y)+F(0)9(-y) + 6f u , (13) 

where 9 is the Heaviside step function and the constants 
F(Q) and F(0) are the homogeneous solutions of Eqs. 
© and HI U| l at y = ±oo. The matrices F have the form 



where 



F (Q) = 
F 3 (Q) = 



F = t 2 {&qF q + a 3 F 3 ), 



lbV(Q) ' 
Vs (k 2 + Q 2 + 2kl) 
lbV{Q) 



(14) 

(15) 
(16) 



and 



V{Q) = (k 2 +Q 2 + 2kl) (k 2 + 2kl) + iki (17) 

The correction 5f u (k,y) obeys the same Eqs. © and 
HlUp without the right-hand side. It has the form 



Sfu = T2O3/3 + T 2 lTo/o + Tl&lfl. 



(18) 



The first term is the singlet component. The second 
term is the triplet component with zero projection of the 



Cooper pair spin on the z axis. This component arises 
even in the case of a homogenous magnetization of the 
ferromagnet and decays in the F film over the short dis- 
tance The last term in Eq. H18p is the triplet compo- 
nent with the spin moment projection ±1. It arises in the 
case of a nonhomogeneous magnetization and decays over 
a long distance of the order £r. The functions fi(k, y) in 
Eq. H18p can be represented as a sum of eigenfunctions 
of Eqs. 10 and (TJJ, i- e -> 

fiiv) =Y. A * iex v(- K i(Q)y)> aty>0, (19) 
l 

My) =^B lie xp( Ki (0)y), at y < 0. (20) 



The inverse decay lengths Ki(Q) are the eigenvalues of 
Eqs. J3J and (fTOf) (without the right-hand side). The 
equation for ki(Q) has the form (/ = 1, 2, 3) 



Q - 2K 



4 (Qki) 



k 2 - 2kf 



k 2 -Q 2 - 2k 2 \ = 0. (21) 



We assume that the exchange length is the shortest 
length in the problem: 



k 2 3> k 2 ,Q 2 , k 2 



(22) 



Then the eigenvalues ki consist of two "short-range" val- 
ues 



K± w (1 =F «sgnw) k h , 
and one "long-range" value 

k l (Q) « y/k 2 + Q* + 2k*. 



(23) 



(24) 



At y < we have the same k; with Q = 0. 

Calculating the corresponding eigenvectors under as- 
sumption H22I) . in the first order over 1/kh we obtain 



a ±~ta 3 ±, A 1z 



2Q . 
T — A 3 ±, 

K± 



A 



3,0L 



0. (25) 



To be more exact, A 3 l ~ (QKL/ik 2 sgnw) Ail and Aol 
is even smaller. The same relations with Q = hold for 
Bn, which yields 



B ± w =F^3± 



(26) 



for nonzero coefficients. 

The next step is to match solutions (|19f) and 1121 III with 
the help of boundary conditions (| 1 1|) and l|12p and to find 
the coefficients An and Bu . This simple but cumbersome 
calculation is presented in the Appendix. In the consid- 
ered limit of a small exchange length [see Eq. J22J] , the 
coefficients A^l ~ -Bil [see Eq. HA.15P ] that describe 
the LRTC are the largest ones. In this limit the function 
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Fq(Q) has a simple form i|A.12|) . Therefore the magni- 
tude of the LRTC at the interface between a domain and 
a domain wall (y = 0) is equal to 



/i(M) = /i(M) = - 



ifssgnu; Q 
Jbkl kq + n ' 



where for brevity we have denoted 

kq = k l (Q), k = K L (0). 



(27) 



(28) 



Now we return to the real space and analyze our re- 
sults. The spatial dependence of the LRTC at y < is 
[seeEq. J2U] 



fUx,y) 



dk 
2^ 



/ L (fc,0) exp(i/cx) exp(K y) 



(29) 



From this formula one can easily find the asymptotic 
behavior of the LRTC /l(x,z/). For large negative y 
(\y\ ^> l/k^), we expand «o and kq with respect to 
(k/ku,) 2 (since the characteristic k in the integral is of 
the order \/%7/\y\) and obtain 



fUx,y) 



Ibkl y V2n\y\ 
exp (-V2k u \y\) exp (-x 2 k^/V2\y\) 



(30) 



This formula shows that the condensate function fi,(x, y) 
decays exponentially with increasing y, but the charac- 
teristic length is rather large (~ fcj 1 ). 

For comparison, we can calculate the short-range com- 
ponent at the SF interface in the case of a homogenous 
magnetization: /o = —ifssgnaj/2~fi,kh [it follows directly 
from Eq. fTT))) ]. Formula Ipto")! shows that at the SF in- 
terface (x = 0) and at distances \y\ < fcj 1 , the func- 
tion /l(x, y) is of the order fs mm(Q, k^/jbk 2 , which is 
smaller than the amplitude of the short-range component 
by the parameter min(Q, k^/kh- Thus, the interface am- 
plitude of the LRTC is smaller, however it decays much 
slower in space. 

In the domain wall (y > 0) the behavior of the function 
fL(x,y) is nearly the same as at y < [Eq. if Q < 

k u . In the opposite limit Q > k u the function /l(x, y) in 
the domain wall decays faster: fL(x,y) oc exp(—Q\y\). 

Having found the condensate function /, we can calcu- 
late the density of states (DOS) in the ferromagnetic re- 
gion. The DOS, normalized to the normal-metallic value, 
is given by the general formula 



v(e) = - Re Tr (f 3 a i 



(31) 



Using the normalization condition g 2 — 1 and the small- 
ness of the condensate function, we can write the correc- 
tion to the DOS due to the proximity effect as 



5v{e) = 



(32) 



(we consider the region in space where only the LRTC is 
essential). This expression is valid at any Q, both zero 
and nonzero. Equations (|30[1 and l)32|l show that the 
LRTC changes the DOS in the ferromagnet at distances 
much larger than the exchange length 



III. TRIPLET COMPONENT IN 
MULTIDOMAIN SF STRUCTURES 

In this section we study the LRTC in a SF structure 
with a multidomain ferromagnetic layer, still assuming 
infinite thickness of the F layer [Fig.^b) with d — > oo]. 
One can distinguish between two possibilities: (a) pos- 
itive chirality, when the magnetization vector M(y) in 
all the domain walls rotates in the same direction (e.g., 
clockwise), and (b) negative chirality, when the vector 
M(y) in neighboring domain walls rotates in the opposite 
directions [e.g., clockwise in the 2n-th domain walls and 
counterclockwise in the (2n + l)-th domain walls]. Wc 
are interested in the LRTC assuming that the exchange 
length £/j is much smaller than the coherence length 
At distances x essentially exceeding the length £/j only 
the LRTC survives in the F layer. 

We assume that the width of the domains with Q = 
is 2ao and the width of the domain walls (Q ^ 0) is 2oq. 
The origin (y = 0) is located in the middle of a domain 
with the constant magnetization. At x 3> £,h only the 
long-range components of the condensate function sur- 
vive in the ferromagnet. The largest long-range compo- 
nent is the LRTC. At the boundary between a domain 
and a domain wall the solution must satisfy boundary 
conditions and l|12|) . Consider first the case of pos- 
itive chirality. The angle a(y) is then an odd function 
of y, which means that fi(y) is also odd — this general 
symmetry can be demonstrated in Eq. 10. Hence the 
solution for the LRTC is 

hiv) = Asinh(K y), -a < y < a , (33) 
h(y) = Bsmh(K Q (y - a - a Q )), a Q < y < a + 2a Q . 

(34) 

Matching these solutions and their derivatives at y = do, 
we find 



B = -A 



sinh 0q 
sinh On 



QF 



tanh ( 



tanh (?o 



, (35) 



where 0q = kqciq and = K a a . The amplitude of the 
LRTC at y — clq is 



/i(ao) 



QF 



KQ coth 0q + Kq coth 0q 



(36) 



We see that /i(ao) turns to zero both at ag — > and 
ao — > 0. These limits mean that the widths of the domain 
walls and domains are assumed to be small in comparison 
with £x while larger than The case an — implies 
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that we have a domain structure with collinear magneti- 
zation orientation. The case ciq = corresponds to a SF 
structure with continuously rotating magnetization (the 
case studied in Ref . l38l) . In both cases, the LRTC does 
not arise. 

The spatial dependence of the LRTC in the domain 
(\y\ < do), corresponding to Eq. (|36[1 . is given by the 
inverse Fourier transformation 

h(x,y)=J-e /iK)-^. (37) 

Interestingly, the function /l(x, y) turns to zero in the 
center of a domain (y = 0). This means that the DOS 
variation due to the LRTC also turns to zero in the do- 
main center. 

Consider now the case of negative chirality, when the 
M vector rotates in the opposite directions in neighbor- 
ing domain walls. In this case the spatial dependence of 
the function fi(y) in domain walls remains the same as 
before, i.e., this function is an odd function with respect 
to the center of a domain wall. However the spatial de- 
pendence of the LRTC in domains changes drastically: it 
becomes an even function with respect to the center of a 
domain. Therefore this dependence is 



h(y) = Ccosh(K y), 
fi(y) = Dsmh(K Q (y 



-a < y < a , (38) 
a o ~ a Q )) , a <y < a + 2a Q . 

(39) 



From boundary conditions and (|12[1 we find the co- 
efficients C and D, and finally 



fi(ao) = 



QF Q 



kq coth 8q + kq tanh 8q 



(40) 



In this case the LRTC disappears only in the limit 
aQ — > because in this limit one again has a domain 
structure with collinear magnetization orientation and 
very narrow domain walls. 

Another type of SF structures, sensitive to the chirality 
of the vector M, was considered in Refs. [2(] and 0. It 
was shown that the sign of the critical Josephson current 
in a multilayered SF structure depends on chirality. 



IV. FINITE THICKNESS OF MULTIDOMAIN F 
LAYER 

In this section we consider a realistic structure with 
a ferromagnetic layer of finite thickness d, see Fig.^b). 
We again have to solve Eq. Q with boundary conditions. 
The first of them, at the SF interface, is Eq. JSJ and the 
other one, at the free surface of the ferromagnetic layer, 
is 



dj_ 

Ox 



= 0. 



(41) 



Similarly to Sec. [HJ we can continue / to the whole x 
axis, reflecting it with respect to x = and periodically 



continuing from the (—d,d) interval. The advantage of 
this trick is that the boundary conditions are included 
into the equation. Similarly to Eq. 10, we obtain 

a£ + w~ 2pJ ~ lkl sgn K^ 3 ' # cosa(y) 



+ h [5-2, f] sma(y) 



.Hi 
lb 



S(x-2dN), (42) 



AT=-c 



which must be solved at all x in the class of even and 2d- 
periodic functions. A periodic function can be expanded 
into the Fourier series: 

/»=W E etknX f(k n ), k n = ^-n, (43) 



f(k n ) = / e- ik ^f(x)dx. 

J-d 

After this transformation, we reproduce Eq. (7J with 
the only difference that the continuous wave vector k is 
substituted by discrete k n . 

As in Sec. [H] the requirement that f(x) is even, is au- 
tomatically satisfied since the equation contains only fc^, 
while the 2<i-periodicity is guaranteed since we consider 
kn defined by Eq. JEJ|. 

The equivalence to the previous equations allows us to 
directly use the results of Sees. |H] and IIIII obtained for 
the infinite d. The rule is very simple: in the case of 
finite d, all the results of Sees. Inland IIIII for the Fourier 
harmonics are valid if we substitute k by k n . The real- 
space function can then be calculated with the help of 
Eq. gSJ- 

For example, we consider the case of multidomain SF 
structure with the F layer of thickness d. For the case of 
positive chirality, instead of Eq. I|37|l inside of the domain 
(\y\ < ao), we obtain 



1 

2d 



smh(n Q y) 
sinh 6q 



(44) 



where /i(oo) is given by Eq. 1361) . 

The formula l|44|l can be drastically simplified in the 
limit when the F film is thin for the long-range compo- 
nent but thick for the short-range one (i.e., kh ^> l/d^S> 
Q, kaj). In this case, the main contribution is given by the 
n = harmonic, since otherwise k in the denominator of 
Eq. (|36|l become very large. Therefore, Eq. (|44|) yields 



Q 



ifssgnu 

2c?7bfc^ J Kq COth^Q + Kq coth 6*o 

sinh(K y) 



sinh(Koeto) 



(45) 



where we have used Eq. (|A.12|I . The x dependence has 
vanished since the F layer is thin. 
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FIG. 2: Correction 5v(y) (due to the proximity effect) to the 
DOS at the free surface of the F layer in the case of posi- 
tive chirality. 5u is multiplied by a dimensionless parameter 
(d7bfc^) , while y is normalized by Q. The curves are plotted 
at several energies e (normalized by DQ 2 ). The width of the 
domains is ao = 5/Q, while the rotation of magnetization in 
the domain walls corresponds to Qciq = it. We assume the 
limit DQ 2 <C A, which means that the width of the domain 
walls is larger than the coherence length. At the boundary 
of applicability of our approximation, Sv(y) falls into experi- 
mentally measurable range— (see text for details). 

The variation of the DOS in space, Sv(y) = 
— Re /£/2| w _ > _ iE ., differs drastically from the case with- 
out the LRTC: it is almost constant across the layer (no x 
dependence) and equal to zero in the middle of domains 
(y = 0). The condensate function fa (y) decays exponen- 
tially (oscillating at the same time) from the boundaries 
between domain walls and domains with characteristic 
length £ T = yjD/2nT (for the DOS, £ e = y/W/e); see 
Fig.[21for illustration. At the same time, the singlet com- 
ponent in the case of a domain structure leads to an even 
dependence of Sv(y) with respect to the middle of a do- 
main. The characteristic length of the short-range com- 
ponents is £h — y/D/h. This case is realized if 9q — > 0. 

The case of negative chirality is treated similarly and 
inside of the domain (\y\ < ao) we obtain 

ft \ 1 ik n *t i xCOshjKny) 

where /i(ao) is given by Eq. (|40|l . In the limit kt ^> 
l/rf> Q,ku, Eq. (O yields 




yQ 

FIG. 3: Correction Sv(y) to the DOS at the free surface of the 
F layer: comparison of positive (denoted by "+") and negative 
(denoted by "— ") chiralities at two widths of domains, ao = 
3/Q and ao = 4/Q. The energy is e — 0.5DQ 2 . The difference 
between the chiralities becomes less pronounced as ao or e 
increases. 

The dimensionless parameter (cfybfc^) 2 , which multi- 
plies Sv(y) in Figs. [21 and |3 is much larger than unity 
under our assumptions. This means that 5v(y) is very 
small. However, at the boundary of our approximation, 
when the parameter is of the order of unity, we expect 
no qualitative differences to appear. At the same time, 
in this case the DOS correction shown in Figs. |3 and |3| 
already falls into experimentally measurable range: in 
Ref. |4J, the resolution of the DOS measurement was of 
the order 0.002 0.003. Experimentally, it is desirable 
to measure 5v at small energies, since the effect grows as 
energy decreases, see Fig. 

The behavior of the LRTC [and the corresponding cor- 
rection to the DOS 5v{y)] inside of the domain walls is 
similar to Eqs. I)44|l and (|46|) with the main difference 
that the decay length in the y direction is determined 
not by kq but by kq. 

The coordinate dependence of the DOS presented in 
Figs. |21 and 13 differs from the corresponding dependence 
of the DOS variation caused by the singlet component. 
If the thickness of the F film d is less than the exchange 
length then the DOS variation, related to the singlet 
component, in a domain wider than does not depend 
on y. In the case of the LRTC the DOS variation is a 
strongly coordinate-dependent quantity since the LRTC 
arises near the domain walls. 



fa{x,y) 



z/gsgnu 
2d lb k 2 h 



kq coth 9q + kq tanh 6*o 
cosh(K y) 



cosh(K ao) 



(47) 



fc„=0 



The resulting corrections to the DOS, 8v{y), in the 
case of positive and negative chiralities are compared in 
Fig. El 



V. CONCLUSIONS 

We have studied the long-range triplet superconduc- 
tivity or LRTC in the SF multidomain structure, where 
magnetic domains are separated by the Neel domain 
walls. The magnetization vector M is supposed to lie 
in the plane of the F film. We show that in this case 
the LRTC arises at the domain walls and decays in do- 
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mains over a large distance £t = \J D /2tvT (we assume 
the diffusive case). The same length characterize the de- 
cay of the LRTC from the SF interface. Although the 
amplitude of the LRTC is less than the amplitude of the 
singlet component in the F film at the interface by k/kh 
times, it decays much more slowly in the F film (here the 
characteristic value of the wave vector k is of the order 
£t or d^ 1 ). Therefore if the thickness of the F film d 
is much larger than the "exchange" length ^ h = y/D/h, 
which characterizes the decay of the singlet component 
in the ferromagnet, only the LRTC / survives at the 
outer surface of the F film in a SF structure. Its spatial 
in-plane dependence in domains f(y) differs drastically 
from the corresponding dependence of the singlet com- 
ponent near the SF interface. If the vector M(y) rotates 
in the Neel walls in the same direction (positive chiral- 
ity), then / turns to zero in the centers of domains. This 
implies that the DOS variation due to the proximity ef- 
fect 6u(e, y) ~ / 2 (e, y) varies inside a domain turning to 
zero in the middle. Thus the measurements of the DOS 
variation at the outer surface of the F film in SF bilay- 
ers allows one to get an information on the nature of the 
condensate in the ferromagnetic films (singlet or triplet). 

The effects of the LRTC are mostly pronounced at the 
boundaries between domains and domain walls. How- 
ever, the correction to the DOS that we find is small 
since we assume the SF interface of low transparency and 
hence a weak proximity effect. At the same time, as we 
increase the interface transparency and reach the limit of 
applicability of our approximation, the correction to the 
DOS falls into experimentally measurable range, while 
we expect that our results are still qualitatively valid. 
Another way to increase the correction to the DOS is to 
measure it at smaller energies. 

The obtained results provide an insight for the Joseph- 
son effect in SFS junctions with multidomain structure 
in the F layer. If the thickness of the F layer d is much 
larger than the short exchange length the Josephson 
coupling between the S layers is due to the LRTC. In this 
case the local critical current density j c {y) is modulated 
in space, reaching its maxima at the domain walls and 
decaying to the centers of domains. 
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APPENDIX: CALCULATING A u AND B it 

To find the coefficients An and Bu , we match solutions 
()19J) and (|20|l with the help of boundary conditions l|ll|) 



and (jT2"j) . This yields 

A 3+ + A 3 _ + F 3 (Q) = B 3+ + B 3 - + F 3 (0), (A.l) 

-A 3+ + A 3 _ + F (Q) = -B 3+ + S 3 _ + F (0), (A.2) 

_^4 3+ + ^U 5 _+A li =B li , (A.3) 

K+ K- 

-k + A 3+ - k^A 3 - = k + B 3+ + k-S 3 _, (A. 4) 

k+A 3+ - K-A 3 - - QB 1L = -k+B 3+ + K-B 3 ^, (A.5) 

Q [A 3+ - A 3 _ + F (Q)} - k q A 1l = KQ B 1L , (A.6) 



where for brevity we have denoted 

kq = kl{Q), k o = k l (0). 
From Eqs. (|A.4jl and l|A.5jl we find 

B 3 + = -A 3+ + -^—Bi L , 

B 3 - = -A 3 _ - -Q-Bxl. 

2k- 

From Eqs. I|A.1|) and 1A.2() we find 

2(A 3+ + A 3 _)+ Q 1 ' 

2(-A 3+ + A 3 _) + 



(A.7) 

(A.8) 
(A.9) 



2 

Q ( 1 





) B 1L = -SF 3 , 


K +/ 




1 ^ 


(A.10) 


) B 1L = -6F , 


K + j 






(All) 



where 8F 3 = F 3 (Q) - F 3 (0) and SF = F (Q) - F (0). If 
Eq. I|22|) is fulfilled, on the order of magnitude we have 
SF 3 cx (Q/k h ) 2 F and SF cx {Q/k h f{{k 2 + 2kl)/k 2 h )F , 
where 



F 



ifssgnuj 
7bk 2 h 



(A.12) 



It follows from Eqs. (|A.10(> and (|A. that the coeffi- 
cients A 3 ± are smaller than B\l by the parameter Q/k^. 
From Eq. I|A.3(I we find 

A 1L ^B 1L . (A.13) 

From Eq. (|A.6|I we find 

Q{A 3+ - A 3 _) - k q A 1l + QF = k B 1l . (A.14) 

The first term here is small, therefore we finally obtain 

QF 



Ail ~ B 



1L 



HQ + Kq 



(A.15) 



These coefficients determine the amplitude of the 
LRTC in the domain (y < 0) and in the region with 
rotating magnetization (y > 0). 
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